DNA bubble dynamics as a quantum Coulomb problem 
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We study the dynamics of denaturation bubbles in double-stranded DNA on the basis of the 
Poland- Scheraga model. We demonstrate that the associated Fokker-Planck equation is equivalent 
to a Coulomb problem. Below the melting temperature the bubble lifetime is associated with the 
continuum of scattering states of the repulsive Coulomb potential, at the melting temperature the 
Coulomb potential vanishes and the underlying first exit dynamics exhibits a long time power law 
tail, above the melting temperature, corresponding to an attractive Coulomb potential, the long time 
dynamics is controlled by the lowest bound state. Correlations and finite size effects are discussed. 

PACS numbers: 05.40.-a,02.50.-r,87.15.-v,87.10.+e 



Introduction. The dynamics of bubble formation in 
double-stranded DNA (dsDNA) is a problem of high cur- 
rent interest in biological and statistical physics. Un- 
der physiological conditions the Watson-Crick double he- 
lix is the equilibrium structure, its stability effected by 
hydrogen-bonding of base-pairs and stacking between 
nearest neighbor pairs of base-pairs By variation of 

temperature or pH-value double-stranded DNA progres- 
sively denatures, yielding regions of single-stranded DNA 
(ssDNA), until the double-strand is fully denatured. This 
is the helix-coil transition associated with the melting 
temperature T m 

Subject to thermal fluctuations dsDNA spontaneously 
unzips and forms flexible single-stranded DNA bubbles 
ranging in size from a few to some hundred broken base- 
pairs, depending on T and salt conditions 0, As- 
suming that the bubble-forming dynamics takes place on 
a slower time scale than the equilibration of the ssDNA 
strands constituting the bubbles, this DNA-breathing 
can be interpreted as a random walk in the ID coordinate 
x, the number of denatured base-pairs. 

On the basis of the Poland-Scheraga model for DNA- 
melting ^(J DNA-breathing has been studied in terms 
of a continuous Fokker-Planck equation [(| , and of a 
discrete master equation and stochastic simulation 0, . 
In the present Letter we show that the Fokker-Planck 
equation for bubble breathing is equivalent to a quan- 
tum Coulomb problem with a repulsive potential above 
T m and an attractive potential below T m . This mapping 
allows us to discuss DNA bubble dynamics in terms of 
the spectrum of a 'hydrogen-like' system and to derive 
several exact results such as the exact scaling of the bub- 
ble survival behavior and the associated correlations. 

Static and dynamic model. The Poland-Scheraga_free 
energy for the bubble statistics has the form 0, \~ 



We here assume a continuum formulation and imply a 
cutoff for x ~ 1. 70 is the free energy barrier for initial 
bubble formation, jx is the free energy for the disso- 
ciation of x base pairs, and cm a: the entropy loss fac- 
tor associated with the formation of a closed polymer 
ring. The free energy density 7 = 71(1 — T/T m ), where 
T m is the melting temperature. From experimental data 
one extracts approximate values for the parameters. In 
units of kT T with reference temperature T r = 37° C, we 
have 70 ~ 10/cT r , ji w AkT r , and c w 2kT r ; the melting 
temperature for standard salt conditions is in the range 
T m w 70 — 100°C, depending on the relative content of 
weaker AT and stronger GC Watson-Crick base-pairs. 

The stochastic bubble dynamics is governed by the 
Langevin equation with Gaussian white noise £(£), 

% = ~ D ^ + ^ (tt)(t)=2DkT5(t), (2) 

where the kinetic coefficient D of dimension (fcT r ) _1 s _1 
sets the overall time scale of the dynamics: [DkT^ 1 ~ 
/is. With dimensionless parameters fi = c/2kT and e = 
(7i/2fcT)(T/T m — 1), and measuring time in units of /is, 
the Fokker-Planck equation corresponding to reads 



dP d / fi 
dt dx \x 



2 dx 2 



(3) 



Note that close to the physiological temperature T r , /i w 
1 and e w 2(T/T m - 1). 

General results. Eliminating the first order term by 
means of the substitution P = e ex x~^ l P, P satisfies 



dP 



l d 2 P 
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(4) 



This is the imaginary time Schrodinger equation for a 
particle with unit mass in the potential V(x) = fj,(fi + 
m l)/2a; 2 — fie/x + e 2 /2, i.e., subject to the centrifugal bar- 
rier l)/a; 2 for an orbital state with angular momen- 
where x > is the bubble size in units of base pairs. turn fj, and Coulomb potential —fie/x. Introducing the 



& = 7o + l x + c In x, 



2 



Hamiltonian H = — (1/2) d 2 /dx 2 + V(x) and expanding 
P on the normalized eigenstates H^f n = E n ^ n , the 
transition probability P(x,x ,t) from initial bubble size 
Xq to a final bubble size x at time t yields in closed form, 

P(x, x , t) = e< x ^ (—) " V e- E »H n (x)* n (x ).(5) 

Here the completeness of ^ n ensures the initial condition 
P(x, xq, 0) = 6(x — xq). Moreover, in order to account for 
the absorbing boundary condition for vanishing bubble 
size we choose 'l'n(O) = 0. We also note that for a finite 
strand of length L, i.e., a maximum bubble size of L, we 
have in addition the absorbing condition ^> n (L) = for 
complete denaturation. Expression (JSJ) is the basis for 
our discussion of DNA-breathing, relating the dynamics 
to the spectrum of eigenstates, i.e., the bound and scat- 
tering states of the corresponding Coulomb problem • 

The transition probability P is controlled by the 
Coulomb spectrum. Below the melting temperature T m 
(e oc (T/T m — 1) < 0), the Coulomb problem is repul- 
sive and the states form a continuum, corresponding to 
a random walk in bubble size terminating in bubble clo- 
sure (x = 0). At the melting temperature (e = 0), the 
Coulomb potential is absent and the continuum of states 
is governed by the centrifugal barrier alone, including the 
limiting case of a regular random walk. Above the melt- 
ing temperature (e > 0), the Coulomb potential is at- 
tractive and can trap an infinity of bound states; at long 
times the lowest bound state dominates the bubble dy- 
namics, corresponding to denaturation of the DNA chain. 
In Fig. ^ we depict the potential in the two cases e ^ 0. 

(i) Long times for T < T m . At long times and 
fixed x and xq, it follows from (jSJ) that the transition 
probability is controlled by the bottom of the energy 
spectrum. Below and at T m the spectrum is continuous 
with lower bound e 2 /2. Setting E k = e 2 /2 + k 2 /2 
in terms of the wave number k and noting that 
^k{x) oc (kx) 1+ ^ for small kx we have P oc exp(— | e | (cc — 
^o))(a;/a;o)" A1 exp(-e 2 t/2) dkeyLv{-k 2 t/2){k 2 xxv) 1+ ^ , 
and consequently by a simple scaling argument the long- 
time expression for the probability distribution 

P(x, x , t) ~ xxl+^e-^-^e-^h-V 2 ^. (6) 

The lifetime of a bubble of initial size Xq created at 
t = follows from Eq. © by calculating the first passage 
time density (FPTD) as time derivative of the survival 
probability, W(t) = - J™ dxdP/dt [13 , or, via Eq. @, 
W(t) = {l/2)[dP/dx+{2fi/x- 2e)P] x=0 . This produces 

W(t) ~ xl+^e^e-^r 3 / 2 -^ (7) 

Below T m , e < and the FPTD W(t) decays exponen- 
tially. The characteristic time scale is set by 

r = 2/e 2 oc (T m - T)~ 2 , (8) 
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FIG. 1: Schematic of the potential V(x) — e 2 /2 of the 
Schrodinger Eq. Q. a) T < T m : The potential is repulsive, 
yielding a continuous spectrum. The bubble fluctuations cor- 
respond to a Brownian walk process in bubble size x before 
collapse at x = 0. b) T < T m . The potential is attractive and 
can trap a series of bound states. At long times the lowest 
bound state indicated in the figure controls the behavior. The 
bubbles increase in size leading to denaturation. 

that diverges as one approaches T m . Finally, from Eq. © 
we infer that P(x,t) oc exp(— ci|e|(x + C2 1 e j i ) ) with con- 
stants Cj > 0, indicating that the profile of the distribu- 
tion has a drift ~ |e| towards bubble closure at x = 0. 

At T m (e = 0) the FPTD falls off like a power law, 
W(t) oc t~ a with scaling exponent 

a = 3/2 + fi. (9) 

The parameter /i = c/2kT (with // w 1 at T w T T ) is as- 
sociated with the entropy loss of a closed polymer loop. 
Ignoring the logarithmic entropic effects (// = 0) we ob- 
tain a = 3/2, characteristic of an unbiased random walk 
0- From Q we also conclude that the mean bubble 
lifetime scales like 

7"mcan 

ct x /\e\ (x x (T m — T) 1 , (10) 

that diverges as the temperature is raised towards T m . 

(ii) Long times for T > T m . Above T m (e > 0) 
the transition probability P is controlled by the low- 
est bound states in the attractive Coulomb potential. 
For the discrete spectrum we have E n = e 2 /2(l — 
(/i/(/i + n)) 2 ), n = 1,2,.... The lowest state is thus 
given by E\ = e 2 ([i + l/2)/(/i + l) 2 , and the corre- 
sponding nodeless normalized bound state by *&i(x) — 
Ax 1+ ^ exp (— nex/(n + 1)) with normalization constant 
,4 2 = [2/xe/(/x + l)] 2 ^+ 3 /T(2^ + 3) 11]. This bound 
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FIG. 2: Lowest bound state of the attractive well in V(x) for 
T = 10Tr (full line) and T = 2.2T r (dashed); with T m = 2T r . 
Inset: log-log plot showing the power-law increase at x = 0. 



state is localized at ~ ^/(T — T m ) and thus recedes to 
infinity as we approach the melting temperature. In 
Fig. [5] we depict the lowest bound state <Fi. From J5J we 
have P(x,x ,t) ~ e e{x ' Xo) (x/x Q y"e~ Elt ^i{x)^i(xo), 
and we note that the dominant contribution to the distri- 
bution originates from the region where the bound state 
peaks, i.e., at ~ 1/(T — T m ). Inserting in (JSJ we obtain 

P(x,x ,t) = A 2 xx 1 +^e^^ + ^ x -^ 1+2 ^ 

Xf ,-f ! (l+2 /J )i/2(l+ fl ) 2 j 

after some reduction. Note from (|llfl that the profile of 
the distribution drifts towards larger bubble sizes with 
velocity ~ e. The associated FPTD becomes 

W(t) = J 4 2 (l/2 + /i)4 +2fl e- £Ko(1+2 ^ /(1+ ^ 



xe 



2 *(l+2^)/2(l+^) 2 



Exact result at T m . At the melting temperature (e = 0) 
the bubble dynamics problem is equivalent to the case 
of a noisy finite-time singularity studied in Ref. [l3| . 
The eigenstates of H are Bessel functions, ^k(x) = 
(kx) 1 / 2 J\/2+fi{kx), and we obtain after inserting in (J5J 
the distribution 

P(x,xo,t) = _ 1/2 _ M / dke~ k t/2 k 2 J 1/2+tl (kx) 

Xq J 

xJi/2 +f i(kx ), (12) 
or, by a well-known identity [l^|. the explicit expression 



P(x,Xo,t) 



XqJ 

^h/2+ t ,{xx /t), 



(13) 



where I v is the Bessel function of imaginary argument 
[Ti |. From we infer the FPTD 



W(t) 



9 1+2/4 



r(i/2 + /x) 

whose maximum at t = ccq/(3 + 2/i) assumes the value 



(14) 



r(l/2 + /x) V3 + 2/i 



2e 



-3/2-jU 



,-1/2 



(15) 



FIG. 3: Bubble lifetime distribution W(t) from Eq. dJ for 
T m = 2T r (full line) and T m = 10T r (dashed). Inset: log-log 
plot of the power-law behavior at long t, with slopes -2 and 
-1.6, as indicated by the straight lines. 



In Fig. we show the zero-size bubble distribution for two 
different melting temperatures corresponding to different 
power-law tails of the FPTD W(t). 

Correlations and finite size effects. In typical experi- 
ments measuring the fluorescence correlation of a tagged 
base-pair, DNA bubble dynamics can be measured on 
the single molecule level [l5j. The correlation function 
C(t) is proportional to the integrated survival probabil- 
ity C{f) oc Jq P(x, xq, t)dx , where L is the chain length 
lia . From the definition of the FPTD we see that indeed 
C(t) = 1 - /J W(t)dt. We find three cases: 

(i) Below T m (e < 0) we obtain from Eq. © C(t) = 
l-xl +2fi e^ Xo J *e- £2t '/ 2 (t')- 3 / 2 -'W, or, in terms of the 
incomplete Gamma function 7 |l4j . 

C(t) = l-4 +2Al el e l a;o (e 2 /2) 1 / 2 +' J 7(-l/2-/i,e 2 t/2). 

(16) 

With 7(0, x) — T(a) — x a e x for i-»cowe obtain 
C(t) = const. + xl+^e^e^t-^ 2 -^-^ 2 (17) 

for large t. Note that the basic time scale of the corre- 
lations is set by e~ 2 oc (T m - T)~ 2 . For t < e~ 2 the 
correlations show a power law behavior, C(t) oc i~ 3 / 2_M ; 
at long times t e~ 2 the correlations fall off exponen- 
tially. The size of the time window showing power law 
behavior increases as T m is approached. In frequency 
space the structure function C(u>) — J exp{iu>t)C[t)dt 
has a Lorentzian line shape for \uj\ <C e 2 , and power law 
tails for \uj\ ;§> e 2 : 



x i+^ e \e\x ( w 2 + ( e 2/2) 2 ) _1 , for |w| « e 2 , 



for > e 2 . 

(18) 

(ii) At T m (e = 0) the exact expression for the FPTD 
(fHlfl combined with relation C(t) — — J Wdt yield 



C(t) = 1 



r(l/2 + y,x 2 /2t) 

r(i/2 + /x) ■ 



(19) 
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For short times t — > the behavior 

lj r(i/2 + A1 ) t 

obtains, while in the long time limit t — > oo, 

2(a;2) 1 /2+M 



C(t) 



(i + 2/*)r(i/2 + m) 



t -l/2-M. 



(20) 



(21) 



(iii) Above T m (e > 0) the DNA chain eventually fully 
denatures, and the correlations diverge in the thermody- 
namic limit. We can, however, at long times estimate 
the size dependence for a chain of length L. From the 
general expression @ we find 



C(t) 



-E,,t, 



V n (x ) / e ex x-^ n (x)dx(22) 



With the lowest bound state ^i(x) = 
Ax 1+fJ - exp (— [iex/(/j, + 1)) of energy E\ = 
( 2 {ji + l/2)/(/i + l) 2 , and after integration over x, 
we obtain 

C{t) OC ^2 e - £ xo(2 AI +l)/( AI +l) e - e 2 (( M +l/2)/( AI +l) 2 )t a ,l+2A I 



x (1 + M )e- 2 [l + (Le/(1 + /x) - l)e ei/(1+ ^] 



(23) 



The correlations decay exponentially with the time con- 
stant ~ £~ 2 (m + l) 2 /(2/i + 1). In frequency space the 
structure function has a Lorentzian lineshape of width 
~ e 2 (2/x + l)/(/it + l) 2 , and for the size dependence one 
obtains 



C(t) 



Le eL/(i+ M ) ; for e L/(l + /x) > 1, 
Le/(l + /i), foreL/(l + /i) « 1 



(24) 



Note that close to T m the correlation function C(t) cx L. 

Summary and conclusion. We demonstrated that the 
breathing dynamics of thermally induced denaturation 
bubbles forming spontaneously in double-stranded DNA 
can be mapped onto the imaginary time Schrodinger 
equation of the quantum Coulomb problem. This map- 
ping allows to calculate the distribution of bubble life- 
times and the associated correlation functions, below, at, 
and above the melting temperature of the DNA helix- 
coil transition. Moreover, at the melting transition, the 
DNA bubble-breathing was revealed to correspond to a 
one-dimensional finite time singularity. 

Our analysis reveals non-trivial scaling of the first pas- 
sage time density quantifying the survival of a bubble 
after its original nucleation. The associated critical ex- 
ponents depend on the parameter /i stemming from the 
entropy loss factor of the flexible bubble, and therefore on 
the ratio T r /T of reference and actual temperature. This 
correction through /i decreases with increasing T. FPTD 



and correlations also depend on the difference T/T m — 1, 
and therefore explicitly on the melting temperature T m 
(and thus the relative content of AT or GC base-pairs). 
We also obtained the critical dependence of the charac- 
teristic time scales of bubble survival and correlations on 
the difference T — T m . The finite size-dependence of the 
correlation function was recovered, as well. 

The mapping of the of DNA-breathing onto the quan- 
tum Coulomb problem provides a new way to investigate 
its physical properties, in particular, in the range above 
the melting transition, T > T m . The detailed study of 
the DNA bubble breathing problem is of particular inter- 
est as the bubble dynamics provides a test case for new 
approaches in small scale statistical mechanical systems 
where the fluctuations of DNA bubbles are accessible on 
the single molecule level in real time. 
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